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ABSTRACT 

We use the photoionisation code CLOUDY to determine both the location and the 
kinematics of the optical forbidden, high ionisation line (hereafter, FHIL) emitting gas 
in the narrow line Seyfert 1 galaxy Ark 564. The results of our models are compared 
with the observed properties of these emission lines to produce a physical model that 
is used to explain both the kinematics and the source of this gas. The main features 
of this model are that the FHIL emitting gas is launched from the putative dusty 
torus and is quickly accelerated to its terminal velocity of a few hundred km s _1 . 
Iron-carrying grains are destroyed during this initial acceleration. This velocity is 
maintained by a balance between radiative forces and gravity in this super-Eddington 
source. Eventually the outflow is slowed at large radii by the gravitational forces of and 
interactions with the host galaxy. In this model, FHIL emission traces the transition 
between the AGN and bulge zones of influence. 

Key words: galaxies: individual: Ark 564, galaxies: Seyfert, line: profiles, radia- 
tion mechanisms: general 



cussion). Further evidence that FHILs arise from an 'inter- 
mediate' region (i.e. between the BLR and NLR) is provided 
by integral field observations of the near infrared FHILs in 
two of the closest AGN : the C irci nus galaxy and NG C 4051 
jMueller Sanchez et al.ll2006l and I Riff el et al.ll2008l , respec- 
tively). These studies showed that these FHILs were being 
emitted at up to 75 pc from the nucleus. 

In their study of the FHILs, iFerguson et all (Il997l ) 
(hereafter F97) used photoionisation models to show that 
a typical AGN spectral energy distribution (hereafter SED) 
could produce these lines in clouds covering a wide range 
of distances from the central engine; from just outside the 
BLR to the extended NLR. However, as their study focussed 
on the relative positions of the emitting regions of various 
FHILs it did not address the issue of the source of these 
clouds and, in turn, explain the kinematics of these clouds 
evident in the profiles of the lines. 

Here, we extend the work of F97 by presenting the 
results of a new set of photoionisation models that, when 
combined with the observed properties of the FHILs, pro- 
duce a self-consistent picture of the sourc e and kinematics of 
the F HIL emitting clouds around AGN. iDorodnitsvn et all 
(2008) recently used hydrodynamic models to show that a 
wind launched from the inner edge of an X-ray illuminated 
torus would reach the velocities observed in the absorption 



1 INTRODUCTION 

The presence of FHILs in the spectra of AGN is a signature 
of the high energy processes occurring in these objects. It 
, has been argued that these lines are a direct response to 
the strong ionising continuum of AGN as they are produced 
by high ionisation states with potentials >100 eV, although 
collisional excitation is still not ruled out as a possible pro- 
duction mechanism (see lPenston et al.ll 19841). 

Early studies of FHILs (e.g. lOsterbrockl Il98ll . 
iPenston et ail Il984 |Pe Robertis fc Osterbrockl 119841 ) ~re- 



vealed that their profiles differed from the more prominent 
low ionisation forbidden lines (such as [O IIIJA5007 and [N 
IIJA6584) and the permitted broad lines (e.g. H/3 and Ha). 
In addition to having widths between these two extremes, 
they were often observed to be blueshifted with respect to 
the low ionisation lines or to possess an asymmetric blue 
wing, suggesting a net flow of the gas producing these lines. 
Irrespective of the interpretation of the source of these kine- 
matics, it is clear that at least part of the FHIL emission 
arises from a region of the AGN distinct from the traditional 
broad line and narrow line regions (BLR and NLR, respec- 
tively; see iMullanev fc W ard 2008 for a more detailed dis- 
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Figure 1. The two S EDs used in our mo dels. The green line 
represents SED-R04 of lRomano et al.l [120041 ',. Points surrounded 
by squares have been extrapolated from the observable portions 
of the SED. The red line represents SED-X constructed using 
optical/UV/X-ray data obtained with the XMM telescope. Dot- 
ted lines indicate the ionisation potentials of the species discussed 
here. Lines showing the 13.6 eV and 1 keV are included for refer- 
ence. 

lines of the so-called warm absorber (<500 km s _1 to 1000 
km s _1 ). As the velocity shifts associated with this feature 
are similar to those of the FHILs we investigate whether 
this same wind could be responsible for the observed shifts 
of these lines. 

Throughout this Letter we use the convention that [Fe 
VII] refers to the [Fe VIIJA6087 line (ionisation potenial: 100 
eV), [O III] to [O III]A5007 (35 eV), [Fe X] to [Fe X]A6374 
(235 eV) and [Fe XI] to [Fe XI]A7892 (262 eV). 



2 THE MODELS 

We base our models on the observed emission lines and 
continuum properties of the NLS1, Ark 5 64, as it is 



known to have very stroin 
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2008 


). Calculations were performed with version 08.00 of 



CLOUDY, last described by iFerland et al.l (|1998T ). to predict 
the emission lines shown in the optical spectrum of Ark 564 
take n with the ISIS spectro meter on the WHT telescope 
(see iMullanev fe WardTl2008l for details regarding these ob- 
servations) . 



2.1 Input SED 

In their stud y of t he low ionisation lines of Ark 564 
iRomano et~ai1 (|2004T ) assembled an SED (hereafter SED- 
R04) that was constrained by their observations between 
the mid-IR-UV and 1-10 keV. They extrapolated the UV 
and soft X-ray continuum to form an SED that peaks at 
^40 eV in vF v space (see fig. [1} . This was used as our ini- 
tial input SED. To determine the sensitivity of our results 
to the shape of the strongly absorbed portion of the SED, 
we assembled our own based on XMM observations of Ark 
564 taken on the 17 th June, 2000 (hereafter SED-X). This 
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Figure 2. [O III], [Fe X] and [Fe XI] profiles decomposed into 
gaussian components, one of which has the same kinematic pa- 
rameters as the fit to the [Fe VII] line. This places upper limits 
on the flux produced with the same kinematics as the [Fe VII] 
line. 



Table 1. Component parameters for the emission lines. The [Fe 
VII] line was fit with a single component. The other line fits were 
then forced to contain a component with the same FWHM and 
velocity shift. 



FHIL 


Flux" 


Shift 6 (km s" 1 ) 


FWHM (km s- 1 ) 


[O HI] 


6.2 


-110 


507 




4.3±0.3 


6±1 


133±2 




8.6±0.3 


-52±2 


193±1 


[Fe VII] 


1.00±0.03 


-110±3 


507±7 


[Fe X] 


1 


-110 


507 




l.liO.l 


-494±13 


620±30 


[Fe XI] 


0.7 


-110 


507 




l.liO.l 


-510±12 


750±28 



"Component flux relative to [Fe VII]. 
^Velocity shift; negative values are blueshifts. 



SED peaks at approximately the same photon energy, but 
has a luminosity at this energy roughly twice that of SED- 
R04. Because the major difference between the two SEDs 
occurs within the unobservable section of the electromag- 
netic continuum, SED-X conforms to the observed fluxes in 
the optical and X-ray wavebands. It should be noted, how- 
ever, that SED-X overproduces flux at 1000 A (3xl0 15 Hz) 
as compared to SED-R04 by a factor of ~2, although this 
could be explained in terms of increased reddening which is 
severe at these wavelengths. We normalised both SEDs to 
hav e a luminosity of 2.64 x 10 43 ergs s _1 at 7000A (4.3xl0 14 
Hz;H 

omano et al.ll2004r ) resulting in total ionising luminosi- 
ties of 4.8xl0 44 ergs s" 1 and 1.1 xlO 45 ergs s" 1 for SED-R04 
and SED-X, respectively. This corresponds to Eddington ra- 
tios of L/L E dd = 1-8 and 3. 9 (respectively) as suming a black 
hole mass of 2.6 xlO 6 M (|Botte et alj|20u4 l 
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Figure 3. Contour plot showing the regions of the parame- 
ter space at which the observed fe-ratios (see i]2,3|l are pro- 
duced (left column) and the [Fc VII] and [Fe X] lines are most 
strongly emitted (right column). The dotted blue contours show 
the log(terminal velocity/km s" 1 ) of the clouds at each point in 
the parameter space. An electronic version of this plot is available 
at www.dur.ac.uk/j.r.mullaney. 



2.2 Cloud geometry, density and composition 

In their work on the line emitting regions of AGN F97 
showed that FHIL emission peaks along ridges of constant 
ionisation paramater, U (= L ion / ' 4mR?nHhuc; where R is 
the distance of the emitting cloud of hydrogen density tih 
from a continuum source with an ionising luminosity of Li 0n ; 
all other symbols have standard meanings). More specifically 
they noted that the optical iron FHILs all reached a max- 
imum luminosity around R ~ 10 18 - 10 19 cm (0.3 - 3 pc) 
and uh ~ 10 5 — 10 7 cm -3 for a source of ionising luminosity 
Lion — 5 x 10 43 ergs s _1 . These values act as a rough guide 
to the location of the FHIL emitters in Ark 564. However, 
as we include radiative driving in our models which tends to 
reduce the density of the emitting clouds and the ionising 
luminosity is greater in Ark 564 we chose to model a large 
portion of the parameter space covering 15 < log(R) < 20.5 
and 4 < log(n,H) < 11 (in steps of 0.5 for both R and uh) 
to ensure we cover all the possible FHIL emitting regions. 

We have modelled our clouds as being dust free. F97 
showed that the inclusion of dust in their models caused the 
FHILs to be suppressed by up to a factor of 25 compared to 
the dust free models. This is largely the result of absorption 
of the ionising radiation by the dust grains as well as the 
depletion of gas phase Fe onto the dust grains. Finally, we 
assumed solar metallicities (our exact ratios are the same as 
those given in F97) 



2.3 Cloud dynamics 

Any model designed to reproduce the observed line ratios 
should also be able to reproduce the velocity structure of 
the FHIL emitting region. To account for this we included 



radiative driving in our models by using the 'w ind' command 
in clo udy, the physics of which is described in lHennev et al.l 
(2005). In allowing the clouds to be radiatively driven, we 
drastically change the physical parameters of the emitting 
gas as compared to models involving only stationary clouds. 
This has corresponding effects on the modelled line ratios. 

By conservation of mass, the density and velocity of a 
continuous column of gas is related such that A-ktqPoVq = 
A-Kr\piv±, where po,i and «o,i represent the density and ve- 
locity at two radii, ro,i. In including radiative driving we 
therefore decrease the density of the emitting clouds as a 
function of depth, allowing ionising radiation to penetrate 
deeper into the cloud. 

The profiles of the emission lines were used to further 
constrain our models. By decomposing the profiles of the 
iron FHILs and [O III] into gaussian components we can 
place upper limits on the luminosity of the lines emitted from 
a region described by a single set of kinematic properties. 
To do this we fit the [Fe VII] line with a single gaussian (a 
second component was found to be statistically unnecessary) 
and force the fits of the other emission lines to contain a 
component with the same velocity shift and FWHM (507 
± 7 km s _1 and -110 ± 3 km s _1 , respectively; see fig. [2] 
and table [1]). From this analysis we note that our models 
of the [Fe VII] emitting gas must not produce [O III] at 
luminosities greater than ~6 times that of [Fe VII] , while the 
[Fe X]/[Fe VII] and [Fe XI] /[Fe VII] ratios of kinematically 
similar components (hereafter referred to as fc-ratios) should 
be no more than unity and 0.61, respectively. To account for 
this and to reduce computation time our model calculations 
were stopped when [O III]/[Fe VII] ^6. Furthermore, we note 
that fig. [2] shows little or no [Fe X] or [Fe XI] flux being 
produced at blueshifts lower than the [Fe VII] component. 

Because photons that are most capable of producing 
Fe 6+ to Fe 10+ are readily absorbed, FHILs will be produced 
nearer the front face of the clouds (i.e. facing the continuum 
source) while the lower ionisation lines will be preferentially 
produced at larger depths. By placing an upper limit on the 
flux of the [O III] emission from the FHIL emitting cloud 
we can therefore limit the size of the modelled cloud. 

Other important stopping conditions were included to 
prevent the clouds from becoming geometrically thick or op- 
tically thick (AR/R < 1 and N H < 10 24 cm -2 ). Following 
F97 the calculation was stopped if the temperature of the 
cloud dropped below 3000K to prevent fragmentation. In 
practice the [O III]:[Fe VII] ratio always became >6 before 
this stopping criterion was met. 



3 RESULTS 

In fig.|3]we present contour plots showing the regions of the 
parameter space at which SED-R04 and SED-X produce the 
observed fe-ratios (defined in W2.'S\i and the [Fe VII] and [Fe 
X] lines are most strongly produced. The [Fe XI] contours lie 
close to the [Fe X] contours and have been removed for clar- 
ity. Also shown on these plots is the velocity the modelled 
clouds reached by the end of the calculation. This could rep- 
resent either the true terminal velocity of the cloud or the 
velocity that the cloud reached at the point when the cal- 
culation was halted, e.g. when the [O III]/[Fe VII] fc-ratio 
reached 6 (these two values may not be the same). Further- 
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Figure 4. Emission line luminosity as a function of gas velocity 
for a cloud at _R = 10 19 cm with njj = 10 6 cm -3 illuminated by 
SED-X (log(C7) = -0.4) 



Figure 5. Same as fig. [4] but for clouds at R = 10 17 cm and 
n H = 3x fO 9 cm" 3 (log(C7) = 0.1) 



more, the FHILs are not necessarily strongest at the termi- 
nal velocity- in general they peak at lower velocities (see figs. 
|4]and[5]). The velocity characteristics of fig. [3] are a result of 
the fact that when the clouds are optically thick essentially 
all the incident photon momentum is transferred to the gas. 
At higher cloud densities the total mass of the clouds will 
increase and hence their terminal velocities decrease. 

Both SEDs are capable of producing the observed fc- 
ratios along a line of constant log((7) = —0.4 stretching from 
log(njj)=5.0, log(i?)=19.5 to log(n H )=10, log(i?)=17.0. 
This suggests that when kinematics are included in the mod- 
els, the precise shape of the illuminating SED has little effect 
on these ratios. However, since both SEDs were produced to 
match the observed emission of a single object, this does 
not rule out the possibility that FHIL ratios will be affected 
by the large differences in overall SED shape between dif- 
ferent AGN. The influence of the SED shape on the relative 
strengths of the FHILs will be the focus of a later paper. 

Neither SED is capable of (radiatively) driving the emit- 
ting clouds to the observed blueshifts of the FHILs in the 
regions of the parameter space where these lines are most 
strongly produced. Even the increased EUV flux of SED-X 
can only produce velocities of ~30 km s _1 and ~100 km 
s _1 for the [Fe VII] and [Fe X] lines (respectively) compared 
to the observed line shifts of 110 km s _1 and 494 km s — 1 . 
As explained in i|2.2l our models do not include dust which 
would lead to increased terminal velocities (see §4). 

In figure |4] we show how the FHIL luminosities vary as 
a function of velocity through an emitting cloud at i?=10 19 
cm and rifr=10 6 cm -3 . This shows that if a cloud were to 
undergo acceleration from 10 km s _1 in this region of the 
parameter space, the [Fe X] line would be overproduced rel- 
ative to the [Fe VII] line at low velocities and the [Fe VII] 
line would be observed at higher blueshifts. This is the op- 
posite of what is seen in the spectra (see fig. [2} and implies 
that the FHIL emitting clouds cannot be undergoing ac- 
celeration in the regions of the parameter space at which 
these lines are produced at the observed fc-ratios. Instead, 
they must have obtained their observed velocities prior to 
reaching these regions. 

Although regions of the parameter space with log([/) = 
—0.4 can explain the low blueshift [Fe X] component and the 
observed fc-ratios, it cannot account for the higher velocity 
[Fe X] component. A possible source of this emission are 
regions of the parameter space with log((7)=0.1 (see fig. [5}. 



Here the luminosity of the [Fe X] and [Fe XI] lines peak 
at velocity v ~160 km s . However, the integrated [Fe X] 
and [Fe XI] flux at v ^110 km s _1 correspond to 35% and 
70% (respectively) of the total line flux. We would therefore 
expect to see significant flux at low blueshifts even at this 
higher value of the ionisation parameter where the clouds 
are most strongly driven. 

In placing limits on the [O III] flux (relative to [Fe VII] 
flux) we are able to derive upper limits on the thickness, 
AR, of the clouds before this line dominates (see i)2.3[) . We 
find that in the case of the [Fe VII] emitting clouds with 
log([/) = —0.4, AR is typically < 1% of R, correspond- 
ing to a column density Nh ~ 10 22 cm -2 . Prior to reach- 
ing this depth, the [Fe VII] emissivity falls by two orders 
of magnitude over the final 70% of the cloud. This implies 
that no FHIL flux is lost in introducing the [O III]/[Fe VII] 
>6.2 stopping criteria. In the case of the clouds described 
in the previous paragraph (with log((7)=0.1) the models 
are stopped when the column density limit (10 24 cm -2 ) is 
reached. As expected, the emissivity of the [Fe X] line drops 
significantly before this depth is reached. 



4 DISCUSSION 

The results from the models described here indicate that 
both SEDs are capable of producing the observed fc-ratios 
along a narrow range of ionisation parameters. This implies 
that the blue wing of the [O III] emission in Ark 564 can be 
produced by the same clouds as are responsible for the [Fe 
VII] line and the low velocity [Fe X] and [Fe XI] components. 
However, these models also indicate that neither SED is ca- 
pable of producing sufficient radiative driving in the regions 
of the parameter space where the FHILs are most strongly 
produced. This suggests that an alternative launching mech- 
anism is required to accelerate the FHIL emitting clouds to 
the velocities seen in these lines' blueshifts. 

In the previous section we showed that the region re- 
sponsible for the observed fc-ratios (log((7) = —0.4) could 
not be undergoing significant acceleration, otherwise we 
would see the effects of this acceleration in the profiles of 
the lines; there would be more [Fe X] flux than [Fe VII] flux 
at low blueshifts. These clouds must be accelerated to the 
observed velocities prior to emitting FHILs. A possible lo- 
cation for this accelerating process is in the previously men- 
tioned [Fe X] and [Fe XI] emitting clouds with log((7)=0.1. 
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However, despite the [Fe X] and [Fe XI] flux peaking at v ~ 
160 km s _1 in these clouds, a significant proportion ol the 
total line flux will, again, be produced at lower velocities 
during the acceleration stage. Furthermore, our models sug- 
gest that even at these increased ionisation parameters, ra- 
diative driving still has difficulty in producing the observed 
line shift velocities (~500 km s" 1 for [Fe X] and [Fe XI]). 

Our results indicate that potential FHIL emitting 
clouds undergo acceleration to the observed velocities prior 
to these lines being produced. Post acceleration, they then 
move outward to larger radii where the conditions are suit- 
able for FHIL productio n (e.g. the torus winds described in 
iDorodnitsvn et alj|200g| ). 

In fig. [3]we show that radiative driving is capable of pro- 
ducing velocities of the order of those seen in the blueshifts 
of the [Fe X] and [Fe XI] lines at 10 17 < R/cm < 10 18 , corre- 
sponding to 0.03 — 0.3 pc. At these radii the energy-density 
temperature, which is roughly the temperature grains would 
assume, is in the range from 800 K to 2500 K for SED-X. 
This straddles the dust sublimation temperature (~1500 K, 
iBarvainisI 19871 ) and is therefore the region at which the inner 
edge of the putative dusty torus would lie. However, some of 
the key FHIL emitting elements (and in particular iron) are 
heavily depleted in dusty regions so the grains must have 
been destroyed by the time the gas becomes highly ionized. 

We therefore suggest that the illuminated face of the 
torus is the origin. We expect dust to have a range of 
grain sizes and a corresponding range of temperatures, with 
smaller dust grains being hotter than larger grains. These 
size and temperature distributions modify the concept of a 
single uniform dust sublimation radius, replacing it with a 
range of radii at which different size dust grains are evapo- 
rated. Further, the AGN luminosity is variable so the posi- 
tion where a given dust temperature occurs changes in re- 
sponse. The large opacity of dusty gas, typically 3-4 dex 
above Thomson, means that the material will be radia- 
tively driven outward while the dust is being sublimated. 
In this scenario we observe the outflowing material (in the 
form of FHILs) whe n the dust extinction becomes small 
|Netzer & Laor][l993l ). The observed coronal line region is 
the result. 

This model facilitates the production of [Fe X] and [Fe 
XI] in clouds across the large range of modelled densities, 
yet showing similar kinematics as the terminal velocity of 
the clouds launched from these radii is largely independent 
of the density of the cloud. 

A launch from the torus in Ark 564 can account for 
the high velocity component of the [Fe X] and [Fe XI] lines. 
Without modification, however, this model would also pro- 
duce the [Fe VII] lines at these same terminal velocities. To 
account for the lower velocity shift of the [Fe VII] line (and 
the corresponding components in [Fe X] and [Fe XI]) a re- 
tardation force must be present at R > 3 x 10 18 — 10 19 cm 
(1 — 3 pc). Beyond this radius our models indicate that more 
than 80% of the [Fe VII] flux, 37% of the [Fe X] flux and 35% 
of the [Fe XI] flux is produced, which is in rough agreement 
with the observed ratios given in table [1] 

We note that the stellar bulge of the host galaxy be- 
comes an important factor at these large radii. An estimate 
of the radius at which the stellar bulge becomes important 
can be determine d by assuming a bulge density distribu- 
tion described by iHernquistl (| 1990T ) . a total bulge mass of 



4.6 x 1O 1O M l|Bian fe Zhaoll2003h and an approximate half- 
light radius for Ark 564 ~1 kpc (determined from 2MASS K- 
band images; this is similar in scale to the ty pical h alf-light 
radii for other SBb galaxy bulges found by iKentl fl985)). 
This reveals that the gravitational potential of the stellar 
bulge becomes important (i.e. greater than the black hole 
potential) at R = 0.03 pc, and the gravitational force of 
the bulge dominates the radiative+gravitational force of 
the black hole/acc retion disk at ~10 pc. More generally 
iDavies et af] (|2007h showed that there is considerable emis- 
sion from young stars within the central few parsecs of an 
AGN-hosting galaxy. It is thus reasonable to consider that 
the kinematics of the outer FHIL emitting region can be af- 
fected by the host galaxy at radii approaching a few parsecs 
from the central engine. 

A number of studies have suggested that winds 
coul d be launched from the inn er e dge of the torus 
(e.g. iKrolik fc Begelmanl Il986l . Il988l and IDorodnitsvn et all 
2008). The results presented here suggest that these same 
winds produce FHIL emission once they reach the optimal 
ionisation parameter for their production. 

GJF is funded by NASA (07-ATFP07-0124) and NSF 
(AST0607028). 
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